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Abstract 

The evolution of a hadronic system after its chemical decomposition is described 
through a model that conserves the hadronic multiplicities to their values at chemical 
freeze-out. The state of the system is found as function of temperature and the cor¬ 
responding baryon density is evaluated. The baryon density at thermal decoupling is 
also computed. 
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1. Introduction 

Thermal approaches have extensively been used to describe the particle multiplicities 
which emerge from high energy collisions [1-16]. The results of such approaches are satis¬ 
factory since they are able to predict a large number of different experimentally measured 
hadronic multiplicities as functions of a few thermodynamic variables, such as temperature, 
volume and chemical potentials. 

The extracted parameters from such approaches are associated with “chemical freeze- 
out” , i.e. the point where the chemical composition of the fireball is fixed. After this stage 
the particles continue to interact until their momentum distribution is fixed, as well. This 
second point is called “thermal freeze-out”. 

Since a set of experimentally measured abundances define the chemical freeze-out point 
these abundances have to remain fixed through the whole cooling process until thermal 
freeze-out. After all, the particles are measured once freeze-out has been completed. 

Various authors have used models of thermally equilibrated relativistic hadronic abun¬ 
dances to determine the chemical freeze-out parameters. In some of these models the hadrons 
are non-interacting particles [2-7] and in others a kind of interaction among them has been 
included [1,8-16]. In this paper the main focus will be on a model of a non-interacting 
hadronic gas formulated in the grand canonical ensemble, called “Ideal Hadron Gas” (IHG) 
[2-5]. It will be argued that in this model, as in the rest of thermal models, it is not pos¬ 
sible to fulfil the requirement that all the particle multiplicities will remain fixed for the 
whole way from chemical decoupling to thermal freeze-out. Considering for example IHG, 
the partition function is expressed as function of (V, T, {A}) where {A} are fugacities associ¬ 
ated with quantum numbers, such as Baryon Number, £>, Charge, Q and Strangeness, S, as 
well as the departure from absolute chemical equilibrium. In this particular case the total 
number of these fugacities is limited to at most five. Not all of them are independent since 
quantities such as <B>,<Q>,<S> have to remain fixed, limiting the total number of 
independent parameters to four. The temperature at thermal freeze-out point is generally 
different from the chemical freeze-out one. So it will be impossible for someone who will 
use the same model at thermal freeze-out point to have all the multiplicities fixed to their 
chemical freeze-out values. The existing free parameters will not be enough. 
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In this work the necessity to have fixed particle numbers will be used to construct a 
thermal model which will determine the evolution of the hadronic system after its chemical 
freeze-out. With the use of this model the construction of the “path” followed by the system 
in the diagram of temperature as function of baryon density will be possible. 


2. The model and its application 

In the context of IHG the grand canonical partition function, formulated in the Boltz¬ 
mann approximation, has the form 


I ' j-fl ■ ■ 

lnZ(V,r,{A}) IHO = VY,*qn,Y. z h„(T) = • « 

i j i j 

where i runs over all hadronic families such as mesons, N Baryons, A Baryons, etc. and 
j represents the specific member of the family with degeneracy factor and mass m %3 . 
XqNi stands for the product of all the fugacities associated with the particular family. These 
fugacities can either be quantum numbers fugacities related to Baryon number, Strangeness, 
etc. or to quark flavour^]. 

Using the above relation one can evaluate particle abundances if one extends the partition 
function by introducing a fugacity \j for every particle. After calculating the particle number 
one has to set X l3 = 1 [17], so again the particle number is only expressed as function of the 
quantum numbers fugacities. 

Now, if someone wishes to keep the particle numbers fixed, e.g. at their chemical freeze- 
out values, it is only natural to use the fugacities X t] , but with the difference that they 
are allowed to be X i3 ^ 1. This model is called Fixed Particle Numbers (FPN) model and 
accordingly the partition function depends on Ay’s 


In Z(V, T,{\}) fpn = VJ2 A h,,Zh,,(T ) 

ij 




H, 


*7 


ij 2vr 2 


9ij m ij 


K, 


,m. 


ij ■ 


( 2 ) 


where A h i:j is product of quantum numbers as well as particle number fugacities^]. The mean 
particle number can be evaluated through the relation 


< N i3 >= A ^ 


d\nZ(V,T,{X}) FPN 


dX i 




( 3 ) 


1 For example, for S - Baryons, A qn would read XbXq 1 Xg 2 "fg or XdX^'y 2 - O ne can l°°k l° r example in 
[16], eq. 14, to find out how the two sets of fugacities are related. 

2 For example, for S(1530)^ Baryon, A h would read A_BAg 1 Ag 2 7 2 AH(i 53 o)-. 
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where {A} 7 ^ X t j means that all fugacities except A^ are considered as constants. 

A thermal model, like IHG, can then be used to extract the set of chemical freeze-out 
parameters {V, T, A b, A s, • ■ •} which best fit a set of experimentally measured multiplicities. 
With these given parameters the numbers < N i3 > of all particles that compose the hadronic 
gas can be calculated. All these numbers have to stay fixed, during the evolution of the system 
after chemical freeze-out, so X tJ will be allowed to become different than one. With the use 
of (3) this requirement takes the form[] 


< N l3 >=< NX VX QNi Z Hij (T) = V'X Hij Z Hij (T') 




A 1/A qn^ZhJT) 

Hi ’ V7.rr l T'\ 


( 4 ) 


In the last equation A u.. may contain quantum numbers fugacities as well the particle number 
fugacity. As it is shown in the Appendix it is not possible to evaluate each quantum number 
fugacity, but this is irrelevant since the full product of fugacities can be calculated. Let me 
point out that all quantum numbers are automatically conserved because they are linear 
combination of the particle numbers. 

The volume of the system at chemical freeze-out, V, on the other hand, need not nec¬ 
essarily stay fixed. If V' at temperature T' is different from V, then all fugacities given by 
(4) depend on a multiplicand factor Ah, which cannot be determined from the constraints 
imposed by the conservation of particle numbers. An additional constraint has to be applied. 
For example conservation of entropy can be assumed^. 

The entropy of the system can be calculated fromf] 


S 


/ a[-rinZ(y,r,Q4)] \ 

^ sr K,m 


( 5 ) 


where /1 represents the chemical potential associated with fugacity A = exp (fJ,/T). Applying 
3 The primed variables in this paper will generally be related to subsequent points of the chemical freeze- 
out point. 

4 A lot of authors assume isentropic evolution of the system, e.g. see [19]. 

5 The symbol of entropy is tilded in order not to be confused with the symbol of Strangeness. K can be 
set equal to one. 
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(5) to the partition function (2)[| the constraint of fixed entropy will read 


^ln Z(V,T, {^}) + VTj2 A 


S = S' 
dZ Hij (T) 


QNi 


V 


d T 




VTY.Iqh^Zh^T) = 


= in z(v, r, {/,'}) + vr y a h ,, az g^ r) - nr y a h „ ^z h>) <r ) . 
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With the use of (4) the last equation becomes 

dz «,A T ) v 
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V'ZhJT) 


( 6 ) 


Setting x = (6) can be solved for x to give[] 


x = exp 


v A 7 I -r ^ \ dZH a (T) w v- a z H ij ( T ) 9Z Hij (T') 

Ziij ^QNi^HijyJ- ) ^\z H .(T')) ' ^ Zij AQNi dT ^ Zij AQN { % h ^p 1 ) dT' 

'Zij ^QNi Z Hii (T) 


(7) 


Equation (7) can be used to evaluate the volume expansion ratio as the system has cooled 
to a temperature T' less than the chemical freeze-out temperature T. With the use of the 
same equation, quantities like the baryon density of the system can be calculated at T'. One 
has to remember that baryon number is also fixed with the imposition of the constraints (4). 
So 


< B' > < B > 


ch 


V <B> 


ch ry~i ch 


n B = 


V' 


V' 


V' 


V 


m 


X 


( 8 ) 


6 The IHG partition function (1), where X l3 = 1, can be used for the evaluation of the entropy at the 
chemical freeze-out point. 

' All hadrons with masses up to 2400 MeV are included in the calculations corresponding to FPN and to 
IHG models. 
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Other constraints which have to be applied to the system and are connected to quantum 


numbers, like 



< B > 
<Q> 


ch 


< |S| >=< \S\ > ch , (9) 


are also satisfied. 

Thus, the contour which is followed by the system after the fixation of its chemical 
composition until its thermal freeze-out can be evaluated. This contour can be defined on 
a ( T,ns ) plane with the use of eqs. (7) and (8) but not on a (T,/ib) plane for the reasons 
explained in the Appendix. 

As an example the model is used to depict the path followed by hadronic systems which 
have been formed at different interactions at SPS after their chemical freeze-out. The chemi¬ 
cal freeze-out parameters used, along with the corresponding references are listed in Table 1. 
From a variety of thermal analyses performed by different authors the particular ones have 
been chosen because they allow for partial strangeness equilibrium (% ^ 1) and they use 
most recent available values for the experimentally measured hadronic multiplicities. The 
values of Table 1 are then used, for each interaction separately, as input to the equations 
< S >= 0 and evaluated through IHG, to determine the rest of the fugacities. 

Thus the whole set of chemical freeze-out parameters (T,/i B ,/ig,/i 5 , y s ) are calculated and 
also the products of fugacities Ain (1) are also set. 

Giving different values to temperature T, equation (8) can be used to calculate the 
corresponding baryon density. The resulting paths for S + S, S + Ag and Pb+Pb interactions 
are shown in Figure 1. For the Pb + Pb interaction the thermal freeze-out temperature is 
calculated in Refs. [21] and [22], For these values baryon density at thermal freeze-out ?z^ er 
can be evaluated. The results are listed in Table 2. The path for Pb + Pb is followed until 
the lower temperature (of the two given in Refs. [21,22]) is reached. 

In order to compare FPN with an IHG model which presents the closest characteristics 
with it, points that are subjected to the constraint that the baryon number and the entropy 
is fixed (< B >=< B > ch and S = S' c/l )[] are also depicted on Figure 1. Let me emphasise 

8 /3 is fixed from the baryon number and charge of the participant nucleons, e.g. see [16]. 

9 The rest of the constraints on quantum numbers, like < S >= 0, etc., are applied as well. 
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that at these points the particle numbers are not conserved. But the IHG points have no 
problem to be depicted on a (T,/ib) plane. This is done in Figure 2. On this Figure there 
is also depicted an “equivalent” value of baryon chemical potential p,s eq for FPN model as 
function of temperature. This chemical potential is calculated through IHG and the only 
connection it has with FPN is that it gives the same baryon density 

1^b(T, flB eq )lHG = nB(T) FPN . (10) 

For comparison with FPN, points that correspond to calculations through IHG for the given 
thermal freeze-out temperatures of Table 2 have also been depicted on Figs 1 and 2. 

Finally in Figure 3 the ratio x = V'/V, where V is the chemical freeze-out volume, is 
plotted as function of temperature for FPN and for IHG (with fixed entropy and baryon 
number). 

3. Conclusion 

After chemical freeze-out the collisions among hadrons that compose the hadronic gas 
can no longer change its chemical composition. Following this requirement an ideal hadron 
gas model (FPN) has been presented that keeps the multiplicity of every particle fixed to the 
value dictated by the chemical freeze-out conditions. In the context of FPN the constraints 
of conservation of quantum numbers are broken up to a larger number of constraints, these 
of conservation of particle numbers. The chemical potentials of quantum numbers are no 
longer “good” variables to describe the evolution of the system. Of course the fugacities 
of particle numbers used as variables in FPN are not “free” parameters. Their values are 
fixed from the given set of the quantum numbers fugacities at chemical freeze-out. So the 
evolution of a hadronic system is described as function of temperature and baryon density 
(after imposing conservation of entropy). This is done for three SPS interactions. 

Following this evolution and using values of thermal freeze-out temperature extracted for 
the Pb + Pb interaction the baryon density at freeze-out is evaluated. As the temperature at 
thermal decoupling for various interactions can be calculated using transverse mass spectra 
or HBT analysis [23] the same procedure can be applied to evaluate the corresponding 
baryon density before free streaming for these interactions. Finally let me point out that the 
necessity to keep the particle multiplicities fixed after chemical freeze-out can be fulfilled to 
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any thermal model, apart from IHG, through the use of the particle fugacities. 

Acknowledgement I would like to thank Professor N. G. Antoniou and Professor 
C. N. Ktorides for reviewing the manuscript and for useful remarks. 

Appendix 

ft will be argued that an ambiguity presents itself when someone tries to evaluate the 
quantum numbers fugacities. Let us suppose that the system is initially in a state which 
is determined by volume V, temperature T, quantum fugacities Xb,Xs, ••• and particle 
number fugacities \j. The question that arises is whether it is possible to determine the 
thermodynamic variables connected to a subsequent temperature T' . Normally the new 
fugacities could be calculated from a set of n equations of the formQ 

<B(V,T,{ A}) >=< B'(V', T 1 , {A}') > , 

< S(V, T, {A}) >=< S'(V', T\ {A}') > , • • • (11) 

which insure for the conservation of the n quantum numbers and a set of m equations of the 
form 

< Nijiy, T, {A}) >=< %(V', T', {A}') > , ■ • • (12) 

which insure for the conservation of the number of the m particle species that are available 
in the hadronic gas. If one tries to solve the above set of the n + m equations one will find 
out that it is impossible to determine all the fugacities. The reason for this is that the n + m 
equations are not linearly independent. When the particle numbers are fixed, automatically 
the quantum numbers are fixed as well (the opposite, of course, is not true). For example the 
equation for the conservation of the baryon number can be expressed as a linear combination 
of equations for the conservation of the number of certain particle species 

< B > - < B' >= (< N n > - < N' n >) + (< N p > - < N' p >) + • • • 

- (< N n > - < N' n >) - (< N p > - < N' p >) - . • • . (13) 

One might think that a way out of the problem of the linear dependency of the equations 

(11) and (12) would be to reduce the number of particle number fugacities by n in which 

10 The primed quantities are connected to T'. 
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case n particle species would be described by only quantum numbers fugacities. Let us 
suppose for simplicity that the hadronic gas is composed only of particles 1 and 2 and the 
only relevant quantum number is B. If one decided to describe the hadronic gas with the 
fugacities A B and A 2 , then two equations would have to be satisfied 


< B >=< B' >V(Z Hl (T)X B + Z H2 (T)\ b X 2 ) = V\Z Hl (T')X' B + Z H2 (T')X' B X' 2 ) (14) 

< N 2 >=< N 2 > V(Z H2 (T)X b X 2 ) = V\Z H2 (T')X' b X' 2 ) . (15) 


Using (15) in (14) one can solve for the final X' B to find 

V_ ZhSJ) VT A' 2 (m 1 /T) 

B V' B Z Ul (T') V'T B K 2 (mi/T') ' 


(16) 


On the other hand if someone had decided to use the set of fugacities X B and Ai he would 
arrive in a similar way to the relation 


„ V , Z h ,(T) VT I< 2 (m 2 /T) 

B V b Z h ,(T') V'T' B K 2 {m 2 /T') ' 1 ' 

It is obvious from comparing (16) and (17)3 ^ ie value of baryon number fugacity 
depends on the choice of which particle number fugacities are kept. This, of course, is 
undesirable. 

Two alternative choices thereby present themselves. The first is to drop the quantum 
numbers fugacities after chemical freeze-out and describe the evolution of the system with 
only the particle number fugacities. The second is to keep the quantum numbers fugacities 
with the ambiguity that accompanies them. In either case the product of fugacities which 
accompany the part of the partition function associated with each particle species has no 
problem to be evaluated. Thus quantities like the baryon density can be calculated. 
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Figure Captions 

Figure 1 Contours (thick lines) that follow hadronic systems after chemical freeze-out on 
(T,ub) plane for 3 interactions at SPS, calculated through FPN (model of Fixed Particle 
Numbers). On the same graph points (dotted lines) calculated through an 1HG model that 
conserve entropy and baryon number are depicted. 

Figure 2 The points of the 1HG model of Figure 1 on the (T, /xg) plane (dotted lines). The 
thick lines represent calculation through 1HG of the baryon chemical potential that leads for 
a given temperature to the same baryon density as FPN. 

Figure 3 The ratio of the volume V' of the hadronic system to its volume V ch at chemical 
freeze-out as function of temperature calculated for the models FPN and 1HG of Figure 1. 

Table Captions 

Table 1 Chemical freeze-out parameters calculated for different interactions at SPS and the 
corresponding references. 

Table 2 Thermal freeze-out temperature calculated in two different references for the Pb ± 
Pb interaction and the corresponding computation of baryon density through FPN. The 
upper errors of baryon density correspond to the upper errors of temperature. The same is 
true for the lower errors. 


Experiment 

T ch (MeV) 

H c £(MeV) 

7? 

Reference 

S+S 200 A ■ GeV 

180.5 ± 10.9 

220.2 ± 18.0 

0.747 ±0.048 

[18,19] 

S+Ag 200 A • GeV 

178.9 ±8.1 

241.5 ± 14.5 

0.711 ±0.063 

[18,19] 

Pb+Pb 158 A ■ GeV 

174.7 ±6.7 

240 ± 14 

0.900 ± 0.049 

[20] 
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Table 1. 


Experiment 

T ther (M eV) Reference 

n^ er ( fm ~ 3 ) 

Pb+Pb 158 A ■ GeV 

Pb+Pb 158 A ■ GeV 

120 ± 12 [21] 

95.8 ± 3.5 [22] 

0 . 09918 : 8 ?! 

0 . 062718:881 


Table 2. 
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